Non-minimal quintessence and phantom with nearly flat potentials 
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We investigate quintessence and phantom dark energy scenarios, in which the scalar fields evolve 
in nearly flat potentials and are non-minimally coupled to gravity. We show that all such models 
converge to a common behavior and we provide the corresponding approximate analytical expres- 
sions for w(fi^) and w(a). We find that non-minimal coupling leads to richer cosmological behavior 
comparing to its minimal counterpart. In addition, comparison with Baryon Acoustic Oscillation 
and latest Supernovae data reveals that agreement can be established more easily and with less 
strict constraints on the model parameters. 
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I. INTRODUCTION 

According to recent cosmological observations, based 
on Supernovae la and cosmic microwave background 
radiation (CMBR) g probes, as well as to WMAP data 
[|[, the universe is experiencing an accelerated expan- 
sion. In order to explain this unexpected behavior, one 
can modify the gravitational theory Q , or construct var- 
ious "field" models of dark energy. The most studied 
models of the literature consider a canonical scalar field 
(quintessence) @, a phantom field, that is a scalar field 
with a negative sign of the kinetic term Q , or the com- 
bination of quintessence and phantom in a unified model 
named quintom [7[. 

In such field dark energy scenarios the potential choice 
plays a central role in the determination of the cosmo- 
logical evolution. However, one can acquire potential- 
independent, general behavior under the assumption of 
slow-roll conditions, which can be embedded in a poten- 
tial chosen to be nearly flat. These potentials have been 
shown to present interesting cosmological features, espe- 
cially in the case where the dark energy equation-of-state 
parameter w is around —1. Although there is not a con- 
crete proof, there are many arguments indicating that in 
the sub-class of "field" dark energy models, nearly flat 
potentials are perhaps a natural and simple way of pre- 
serving w w — 1, either in quintessence Q or phantom @ 
case (although a sufficiently large Hubble friction could 
also lead to such a w for arbitrary potentials). Finally, 
we mention that nearly flat potentials, which keep the 
variations of the scalar fields from their initial to their 
present values small, can also be efficient to avoid un- 
known quantum gravity effects [Toj . 

In [ll| it was shown that in all quintessence models 
with nearly flat potentials, the evolution of w converges 
to a common behavior and can be described by a unique 
expression, which can be approximately given analyti- 
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cally in terms of its present value together with the value 
of dark-energy density parameter at present. In [12| this 
result was extended to phantom models, while the gener- 
alization to quintom scenario was performed in (l3l | where 
a universal expression for w was also extracted, allowing 
for a crossing of —1. Recently, a similar study with a 
minimally coupled, tachyon-type scalar field was imple- 
mented in(I3|. resulting exactly at the same equation of 
state with the one obtained for a standard minimally 
coupled scalar field [III [l2| • 

On the other hand, dark energy models where the 
fields are non-minimally coupled to gravity [la . flq have 
been shown to present significant cosmological features 
[ItJ ■ In our recent work we investigated non- minimal 
quintessence with nearly-flat potential, in the context of 
Brans-Dicke framework. Under a conformal transforma- 
tion it becomes a coupled quintessence model and a uni- 
versal expression for w can be extracted, depending ad- 
ditionally on the coupling parameter. 

In the present work we are interested in studying both 
quintessence and phantom scenarios with nearly-flat po- 
tential, in the general non-minimal framework. The plan 
of the work is as follows: In section [II] we present the 
non-minimally quintessence and phantom models and we 
extract the approximated general solutions for the sub- 
class of nearly flat potentials. In section Hill we compare 
our formulae with the exact cosmological evolution and 
we provide the observational constraints on the parame- 
ters of the model. Finally, section IIVI is devoted to the 
summary of the obtained results. 



II. EVOLUTION OF NON-MINIMALLY 
QUINTESSENCE AND PHANTOM WITH 
NEARLY FLAT POTENTIALS 

Let us construct quintessence and phantom models 
with the field being non-minimally coupled to gravity. 
In order to incorporate both scenarios in a general and 
unified way, in the following we introduce the usual e- 
parameter, acquiring the value +1 for quintessence, and 
— 1 for the phantom case. Throughout the work we con- 
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sider a flat Robertson- Walker metric: 
ds 2 = dt 2 - a 2 (£)dx 2 , 



(1) 



with a the scale factor. 

The action of a universe constituted of a non- minimally 
coupled field <p (canonical or phantom) is 16]: 



(2) 

where k 2 = 8irG, £ is the non-minimal coupling param- 
eter, R = 6(^ + is the Ricci scalar and a dot denotes 
differentiation with respect to cosmological time. The 
Fricdmann equations read: 



:H 2 



94, + Pm, 



(p<p + p<p) + p m (l + W m ) 



(3) 



(4) 



In these expressions, p m is the matter energy density, 
with w m the matter equation-of-state parameter, p^ and 
Ptf, are respectively the energy density and pressure of the 
non-minimally coupled scalar field, given by [lSl [l6|: 



P4, = e\^ 2 + V(<t>) + e&^HU + e3£#V 



(5) 



-e2f(l - 6£)i^ 2 - e3£(l - 8£)i^ 2 + 2^V'(0).(6) 

In such a scenario, dark energy is attributed to the scalar 
field, and thus its equation-of-state parameter reads: 



W4, 



P<j> 



(7) 



where we can equivalently consider the barotropic vari- 
able 70 = Wcj, + I. Finally, the equations close by con- 
sidering the evolution equation for the scalar field, which 
takes the form [Tj| [7j| : 



4> + ZH4> + 6£(iJ + 2H 2 )<j) + EV\<t>) = 0. 
It proves convenient to introduce 

3H 2 

Ptot = P<P + Pm — K~ 



Ptot = V<i> + W m p m , 



and thus to define 



W t ot 



Ptot 



n 2 P4, 



Ptot 3H 2 



(8) 

(9) 
(10) 

(11) 



where O m is the matter density parameter. Therefore, 
the second Friedmann equation can be rewritten as 

k 2 3 
H = - y(ptot +Ptot) = ~2 H i l + Wt °t)' ( 12 ) 
The variable Wtot is an intermediate one and will be useful 
in order to eliminate H (using relation (|12p ) from all 
equations, since, contrary to the minimal coupling case, 
now we face the presence of H in the definition of the 
dark energy equation-of-state parameter (0. 

At this stage, we desire to transform the aforemen- 
tioned cosmological system into an autonomous form 
[l9| . This will be achieved by introducing the auxiliary 
variables: 



K0 



y = 



71' 



and defining A = -^yjjjjj- 

Using these variables, relation (0 gives: 



Sit ee ^-^r = ex 2 + y 2 + eYl£zx + e6^ 5 



(13) 



3H 2 



(14) 



and thus 



Q m = 1 - fi^ = 1 - (ex 2 + y 2 + sl2£zx + e6£z 2 ), (15) 

where fi^ is the scalar field (that is dark energy) density 
parameter. Similarly, substituting (|12[) into ^ expressed 
in terms of (x, y, z) we obtain: 



3iJ 2 



r(l - 4^)a; 2 - y 2 (l + 2\/6£Az) + e4£ 



xz - 



+£l2£V+£6£(l-6£)zVot, (16) 

and then inserting the above formulae into (jlip we can 
explicitly express Wtot, namely 



Wtot = 



e(l - 4g)x 2 - y 2 (l + 2y / 6CAz) + eAjxz + el2gV + w m Q. r , 
1 - e6£(l - 6£)z 2 



(17) 

Inserting this expression for u> tot into (|12p we can elim- 
inate if from (JSj) , (J6j) and thus obtain the dark-energy 
equation-of-state parameter from ([7]) as: 



3 



= l+ws = 2 



x 2 [e- 



e(2 + e3z 2 )g + 18z 2 g 2 ] + 4£xz[2e + 9z 2 £(6g - 1)] + jz{-y/&y 2 \ + 18^ 3 (6g - 1) + e3z[l + 2g + y 2 (6g - 1)]} 



(x 2 + ey 2 + 12£xz + 6z 2 £)[e + 6z 2 £(6(, - 1) 



(18) 



where for simplicity we have considered the dust case 
w m — 0, although this is not necessary. 

Finally, noting that for any quantity F we have F 

tlx cosmological system itself takes the following 



dF 
din 



H 

autonomous form (see also |20( for details): 



dx 
din a 



-3x- 12£z + e 



2 {X - 



6^)(1 



3 3 

+ -y{l + w tot ) 



-76A 2 (r-l): 



Wtot) 

(19) 

(20) 
(21) 
(22) 



v"v 



and w to t is determined by (jTTJ) with 



dy 
din a 

drj 

din a 

dA 
din a 

where T 
w m = 0. 

The autonomous equations (fll?)) - (|2"2"|) correspond to 
the exact cosmological evolution for non-minimally cou- 
pled quintessence or phantom models, with the dark- 
energy equation-of-state parameter wj, and the dark- 
energy density parameter given by ([T5jl and ([13)) re- 
spectively. We mention that when £ = 0, that is in the 
minimal-coupling case, all equations coincide with those 
of PHI. 

In order to acquire analytic solutions and in par- 
ticular to obtain an expression for the experimentally 
accessible quan tity w^(O^), we follow the strategy of 
[H El El Hi. That is > we first differentiate |QSJ| and 
([TJJ with respect to In a, and making use of (119p - (l2"2"l) 

we obtain f'?*' and f . * as functions of x, y. z, A. Then 

d In a a In a ^ ' ' 

using (|18[) and (| 14[) we eliminate the auxiliary variables 
in favor of 7^ and acquiring ,?f and as func- 
tions of 70,^0. Finally, we obtain the exact differential 
equation for ^3^- as 



d~/<t> 
d In a 

d In a 



(23) 



However, comparing to the minimally coupled scenario 
there is a difference. In particular, in the former case 
[ill E2I EH 70 and fi^ depend only on a; and y, and 
thus one can easily acquire the expressions #(70, and 
2/(7^,^0) in order to eliminate the auxiliary variables. 
On the contrary, in the non-minimal coupling case, where 
the richer dynamics is reflected in the appearance of an 
additional degree of freedom, this procedure is not pos- 
sible since one has to eliminate three variables {x, y, z) 



in terms of two (7^ and £1$). Therefore, we make the 
following assumption: Since z in all relations appears al- 
ways multiplied by £, which is in general a small quantity, 
its time variation will not have significant effects. More- 
over, we are considering potentials to be sufficiently flat 
satisfying slow-roll conditions (see below). Under this 
assumption, the variations of the scalar field with time 
will always be sufficiently small. Thus we replace z by 
its average value zq. Doing so, the elimination of x and 
y in terms of 7^ and can be performed, leading to 

the differential equation — /(70, fi^, A), which is a 
rather complicated expression. 

At this stage, similarly to [ll], El EH , we proceed to 
the usual assumptions, based on the nearly flatness of 
the potentials, quantitatively expressed as 



1 dV 
1 d 2 V 

« 1. 



V dtp 



(24) 
(25) 



First, from the definition A 



3 and usin e 

we conclude that A can approximately considered to be 
a very small constant, so that 



A~A - 1 dV 



<j>=cj>a 



«1, 



(26) 



where Ao and <po are the initial values of A and <fi respec- 
tively. This approximation is also consistent with the 
A-evolution equation (f22l) . since d ^ is proportional to 
A 2 , i.e it is very small, justifying that A remains almost 
constant. The second assumption is that, as mentioned 
in the introduction, w$ remains close to —1 throughout 
the evolution, in agreement with observations. Thus, 
I70I <C 1, too. Finally, in the non- minimal coupling case 
at hand we make an additional assumption that was not 
present in minimal coupling models [ll], Ell EH • Namely, 
since the non-minimal coupling constant £ is usually a 
small number, with the conformal value £ = 1/6 being 
the one that accepts a reasonable theoretical justification 
[3ElIE3]; we can assume that ^-Cl, too. 

Under these approximations the complicated differen- 
tial equation — /(t^Q^, A) acquires the following 
simple form: 



27^ 



2 



+ 



3(i-n*)n, 3 (i-n*)v/n^ 



+ e4? 



V s ' 2 7<t> 



3/2 



.(27) 
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Note that in the phantom case 7^ < and e = — 1 while 
in quintessence 7^ > 0, e = 1, thus the combination 
£70 is always positive and the aforementioned equation 
is well-defined. 

Differential equation (|27|) describes approximately the 



cosmological behavior of the model at hand. Note that in 
the limit ^ — ^ it coincides with the corresponding one 
of [111, [13] • It can be transformed into a linear differential 
equation under the transformation s 2 = £7^ and can be 
solved exactly. Thus, the general solution reads: 



7^,(^0) = 1 + = 
1 -a 



1 

= £ 9 







6 6>/2aibCB(n*;i- y.-l + yj + [V3A (1 



20 - e6a/2.zo£] B 



O 3 
2 



3 



- +f)l ■<-*> 



where -B(u; a, 6) is the incomplete beta function defined 
as 



B(u; a, b) 



t a -\i-ty- l dt. 



(29) 



It can be clearly seen that in the limit £ — ► solution 
(I28p coincides with the corresponding expressions of (Til . 
Il2j (for quintessence and phantom respectively), namely: 



= £ 



\2 
A 



1 



n 



__i ±-1 i„ L _ 



4> 



(30) 



In order to express our result in a form more suitable 
for comparison with observations, and similarly to (ill . 



1 12 . Il3l . |18| | , we will use the relation 



l+ffi^-lla 



(31) 



which is the zeroth order solution of ^ = gi(f2^,7^), 
namely 



din a 



3^(1 - fi*) + 0(£ 2 ) + £>(7|) + 0(Ag), (32) 



where fi^o is the value of f2^ at present (ao = 1). We 
mention that (|3Tj) holds only for |7^| <§; 1, that is for 
ss — 1. Finally, note that one could equivalently use 
the redshift z instead of the scale factor a, using a = 
(l + z)'\ 

Using (f3"Tj) solution (f2"g| writes: 



1 + ra^(a) = e- 



1 / [l + - l) a_3 ] (1 - «0o) ' 



e&j2z<£B( [l+(0^ -1) 



1 + (a 3 - 1)^0 j 

+[V3A (1 - 20 - e&V2z £] B f[l + (fi^ 1 - l) a' 3 ] _1 ; 



'2 3 ' 3 1 



1 3 4£ 



3 



(33) 



Relations (f2"5| and (j3"3")) are the main results of the 
present work. They provide w^(O^) and w,p(a), with f^o, 
Ao, zo an d of course £ as parameters, that is they give 
the unique, convergent behavior for the whole sub-class 
of models where non-minimally coupled quintessence and 
phantom fields evolve in nearly flat potentials. In prin- 
ciple, one can eliminate either Ao or zq in terms of w^o 
(the present value of Wd>). However, contrary to the min- 
imally coupling case [lj [l^] this straightforward proce- 
dure leads to a complicated expression for w$(a), and 
thus we prefer to keep (|33p . Finally, we mention that in 
the limit £ — > ([3"3l coincides with the corresponding 
relations of 



III. COSMOLOGICAL EVOLUTION AND 
OBSERVATIONAL CONSTRAINTS 

Having obtained approximate analytical solutions, we 
wish to explore them, examining their cosmological im- 
plications, and to compare them to the numerically- 
elaborated exact evolution for a few different models. 
In figfT] we present the behavior of w,j,(Q ( j > ) for the 
quintessence scenario, in the conformal case, i.e with £ 
being equal to 1/6. The solid curve depicts the w(Q, < f > )- 
behavior as it is given by (|2"5)) with e = +1, while the 
dashed and dotted curves correspond to the exact evolu- 
tion for the potentials V{<f) = (j> 2 and V(4>) — <t>~ 2 respec- 
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FIG. 1: The behavior afw#(Q<j>) for the quintessence scenario 
in the conformal non-minimal coupling case, i.e with £ = 1/6. 
The two curve-sets correspond to Ao — 1 (upper set) and Ao = 
0.5 (lower set). For each set, the dashed and dotted curves 
depict the exact numerical evolution for potentials V(cj>) = (j> 2 
and V(cj>) = 4>~ 2 respectively, while the solid curve provides 
the w^(^l^) -behavior as it arises from our analytical result 
$28\) with e = +1. In all cases the present value of fi^o has 
been fitted to be equal to 0.7 and the initial -value has been 
forced to be very close to — 1 (which yields a relative small 
number for zo, of the order of 10~ ). 



tively. Finally, we present the aforementioned plots for 
two Ao-values, namely Ao = 1 (upper set) and Ao = 0.5 
(lower set). In all cases the present value of fi^o has been 
fitted to be equal to 0.7 and the initial value of has 
been forced to be very close to —1, which yields a rela- 
tive small number for zq (of the order of 10 -5 ) and this 
fact indeed justifies our assumption to replace z by its 
average zq and neglect terms of O(zq). 

As expected, agreement is not satisfactory for A = 1, 
but it improves significantly for smaller values. Further- 
more, since we have neglected terms of order (1 + w^) 2 , 
the resulting deviation between the exact evolution and 
the analytical expression is larger for larger difference 
between u>0 and — 1. These features are valid for com- 
pletely different potentials. We mention that the corre- 
sponding deviations are larger comparing to the minimal 
coupling case of [III > due to the additional errors brought 
in by the C(^ 2 )-approximation, but this difference is rel- 
atively weak due to the smallness of the conformal value 
(£ = 1/6). Finally, in fig. [5] we present the behavior of 
u>0 (a) for all the cases of fig. Q] as it arises from (|33j) with 
e = +1. Again we see that our analytical expression is 
accurate as long as w^a) is not far from —1 and Ao is 
small. 



FIG. 2: The behavior of w<f,(a) for the quintessence scenario 
in the conformal non-minimal coupling case (£ — 1/6). The 
two curve-sets correspond to Ao = 1 (upper set) and Ao = 0.5 
(lower set). For each set, the dashed and dotted curves depict 
the exact numerical evolution for potentials V((f>) = 4> 2 an d 
V(4>) = <j)~ 2 respectively, while the solid curve provides the 
w(a) -behavior as it arises from our analytical result i33t with 
£ = +1. 



From these two figures we observe that when the po- 
tential is nearly flat, then all the examined models con- 
verge to a common behavior (our analytical solutions (|28|) 
and (|33p V Although non- minimal coupling brings an ad- 
ditional degree of freedom, i.e an additional parameter, 
the small value of the coupling (for instance £ = 1/6 in 
the conformal case) downgrades significantly the corre- 
sponding complexity. 

In fig. [3] we show w^, (£7^ ) for the phantom scenario 
with £ = 1/6. In this case, and due to the "inverse" be- 
havior of phantom fields in potential slopes, we expect A 

(defined as A = — ^ ) and thus A too, to be negative. 
This was also incorporated in the minimal coupling case 
[T3 | where the authors defined A straightaway with an 
opposite sign. Furthermore, we mention that the chosen 
potentials are not too stiff, in order to avoid a Big Rip 
which is common in phantom scenario HHH. As wc sec, 
the agreement of our analytical expression with the ex- 
act evolution is better for smaller | Ao | and 1 1 + | . The 
same features are observed in the corresponding u>^(a)- 
behavior shown in fig. [4] 

Having determined the accuracy region of our analyt- 
ical expressions, we compare our result for w,p(a) (given 
by ([33]) ) to supernova Type la and Baryon Acoustic Os- 
cillations data. In this expression for u>^(a) there are 
four parameters, namely f^o> zq and Ao- As we have 
mentioned, we have assumed the conformal coupled case 



0.1 0.2 0.3 0.4 0.5 0.6 
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FIG. 3: (fl^) for the phantom scenario in the conformal 
non-minimal coupling case = 1/6)- The two curve-sets 
correspond to Ao = —0.5 (upper set) and Ao = —1 (lower set). 
For each set, the dashed and dotted curves depict the exact 
numerical evolution for potentials V((f>) = 4> 2 and V(4>) — 
<j>~ 2 respectively, while the solid curve provides the w < f,(Q,^,)- 
behamor as it arises from our analytical result H28V with e = 
-1. 



FIG. 4: w<f,(a) for the phantom scenario in the conformal 
non-minimal coupling case = 1/6). The two curve-sets 
correspond to Ao = —0.5 (upper set) and Ao = —1 (lower set). 
For each set, the dashed and dotted curves depict the exact 
numerical evolution for potentials V((f>) = 4> 2 andV((j>) = (j>~ 2 
respectively, while the solid curve provides the w$ (a) -behavior 
as it arises from our analytical result f33p with s = — 1. 



(£ = 1/6) and we have considered z to be of the order 
of 10~ 5 for our present investigation. While taking dif- 
ferent values for £ will definitely change our result, our 
conclusions below are not so sensitive on zq as long as 
it is sufficiently small. This is also one crucial require- 
ment for our analytical result to be valid. With two re- 
maining parameters Q^q and Ao, one can relate Ao to the 
present day value of the equation of state w^o. Hence 
our aim is to put constraints on the model parameters 
u>0o and f^o- At present, we have the Union08 compi- 
lation of the Snla data which contains around 307 data 
points [2^ . This is world's published first heterogeneous 
SN data set containing large sample of data from SNLS, 
Essence survey, high redshift supernova data from Hub- 
ble Space telescope, as well as several small data sets. 
We use this data set together with the BAO data from 
SDSS (Sloan Digital Sky Survey) [Hj]. The correspond- 
ing confidence contours are depicted in fig. [5] 

We should mention that while fitting our model with 
the observational data, we have to do it separately for 
quintessence and phantom case, since the expressions for 
(1 + w^) and —(l + w$) as a function of scale factor are 
not the same. This is unlike the minimally coupled case 
HHG3, where they are exactly identical. In this regard, 
and from the fitting procedure point of view, one does 
not expect in general to have the confidence contours for 
quintessence and phantom case to match exactly to give 



a complete ellipse, simply because we are fitting two dif- 
ferent w<p(a). On the other hand, since we have switched 
the sign of the kinetic term of the scalar field in order to 
obtain the phantom model, we theoretically expect con- 
tinuity for the exact solution and the contour to cross 
smoothly to the phantom region. Thus, the fact that the 
contours match exactly, as shown in fig. [5j indicates that 
our approximate method works well for both quintessence 
and phantom cases. 

As we observe, and similarly to the minimally-coupled 
cases (ll|, [HI, 03], in order to acquire tight constraints 
on w^o one has to impose strong bounds on fl^o, arising 
from the corresponding constraints on Q m o from struc- 
ture formation data like the shape parameter for the mat- 
ter power spectrum. However, with the present data, 
and although the behavior is richer than the minimally 
coupling scenarios, one cannot distinguish between the 
present model and a ACDM one. 



IV. CONCLUSIONS 

In the present work we have studied dark energy sce- 
narios in which the scalar field is non-minimally coupled 
to gravity and it evolves in nearly flat potentials in order 
for the slow-roll conditions to be satisfied. For complete- 
ness, we have investigated both quintessence and phan- 
tom cases in a unified way, having in mind that the later 
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FIG. 5: la and 2a confidence contours in the plane defined by 
the present-day values of O,^ and w§ (denoted by Q^a and w^o 
respectively). The dashed curves correspond to quintessence 
scenario while the solid curves to the phantom one. Note that 
the contours match exactly at the w^o = — 1 line, although the 
corresponding fitting-expressions are not "symmetric" between 
the two cases. 



paradigm could have problems at the quantum level (the 
relevant discussion is still open in the literature since 
for instance in 24] the authors reveal the causality and 
stability problems while in [2oT ] the authors construct a 
phantom theory consistent with the basic requirements 
of quantum field theory with the phantom fields arising 
as an effective description). 

We have shown that all such models converge to a com- 
mon behavior and we have extracted the corresponding 
approximate analytical expressions for io^(O^) (relation 
(|28p) and for w<p(a) (relation (|33|) ). The smallness of the 
non-minimal coupling parameter (£ = 1/6 in the confor- 
mal case) enhances the validity of our formulae, since it 
downgrades the role of the additional degree of freedom 
comparing to the minimal-coupling case. In particular, 
the errors are typically within 2%, and are smaller for w<f, 
being closer to —1 and for more flat potentials. Moreover, 
a comparison with observations using Baryon Acoustic 



Oscillation and latest Supernovae data reveals that the 
situation is similar to that of a minimally coupled scalar 
field. That is, under the present data the model can- 
not be distinguished from a ACDM scenario. However, 
one should note that in our investigation we have fixed 
the parameter £ to its conformal value 1/6. Keeping £ 
arbitrary will make the parameter space broader, which 
in turn will lead further to less stringent constraints. In 
other words, switching on non-minimal coupling it be- 
comes easier to fit the observational data, at least at the 
background cosmological level. 

Non-minimal coupling leads to richer cosmological be- 
havior comparing to the minimal scenario. For instance, 
although in minimal coupling with completely flat po- 
tentials one obtains the cosmological constant universe 
(i.e w = — 1 = const.) [ll|, [l2j], in the model at hand, 
even in this case one obtains a varying w resulting purely 
from non-minimal coupling. In addition, it is interesting 
to see that non-minimal coupling introduces an "asym- 
metry" between quintessence and phantom cases, which 
is visible straightaway in the energy density and pres- 
sure definitions ([5]) and ([6]). However, even if the expres- 
sions are non-symmetric in the two sides of the phantom 
divide, the observation comparison reveals a non-trivial 
match of the confidence contours at the w^q = — 1 line. 
This feature provides a physical self-consistency test of 
the obtained analytical expressions of the present work. 

Finally, we mention that we have only considered 
observational data related to distance measurements, 
depending only on the background homogeneous and 
isotropic universe. An additional consideration of the 
growth of matter perturbations in these models, can have 
interesting effects and may help to distinguish between 
the present model and the cosmological-constant one. 
This will be our aim in a future investigation. 
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